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We study the performance of a quantum wire spin filter that is based on the Rashba spin-orbit 
interaction in the presence of the electron-electron interaction. The finite length wire is attached to 
two semi-infinite nonmagnetic leads. Analyzing the spin polarization of the linear conductance at 
zero temperature, we show that spin-filtering is possible by adequate tuning of the system parameters 
first considering noninteracting electrons. Next, the functional renormalization group method is used 
to capture correlation effects induced by the Coulomb interaction. For short wires we show that the 
energy regime in which spin polarization is found is strongly affected by the Coulomb interaction. 
For long wires we find the power-law suppression of the total conductance on low energy scales 
typical for inhomogeneous Luttinger liquids while the degree of spin polarization stays constant. 
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I. INTRODUCTION 

The birth of spintronics can be dated back to the dis- 
covery of the giant magnetoresistance effect in 1988,-1*^ 
Since then, many theoretical and experimental studies on 
spin-dependent electronic transport have been performed 
in order to achieve a clear understanding of the under- 
lying physics and to investigate the possibility of fabri- 
cating spintronic devices i^i^i^ One manifest realization of 
such a device is a quantum wire with parameters which 
can be tuned by a set of gate electrodes such that the 
transport of electrons with a certain spin direction is fa- 
vored. The original proposal by Datta and Das^ is to use 
the spin precession in a narrow-gap semiconductor wire 
with spin-orbit coupling between two magnetized leads to 
modulate the current. However, Streda and Seba came 
up with an idea of obtaining spin polarization in a (al- 
most) nonmagnetic system. They considered transport 
through an infinite quasi one-dimensional (ID) quantum 
wire with Rashba spin-orbit interaction (SOI), a mag- 
netic field in the direction of the propagation as well as a 
potential step.^ In their setup, leads with negligible SOI 
and vanishing magnetic field as present in any experi- 
mental realization were not taken into account. 

The goal of the present work is twofold. We first in- 
vestigate, under which conditions spin-polarized currents 
through a finite length ID wire with SOI and parallel 
magnetic field can be achieved, if the coupling to two 
semi-infinite leads is included. We study the influence of 
a potential step and a localized impurity. Such inhomo- 
geneities play an important role for the spin polarization. 

In the Datta-Das setup as well as in the system sug- 
gested by Streda and Seba the spin dependent transport 
properties heavily rely on the strictly ID nature of the 
quantum wire with only one partially filled subband. Al- 
ready in the presence of a few filled subbands the con- 
trol over the spin ceases.— In a few hundred nanometers 
long strictly ID wires the Coulomb repulsion of electrons 
is expected to have a dramatic effect on the physics. 



Such systems cannot be described by Fermi liquid the- 
ory. Instead, the low-energy properties arc captured by 
the Luttinger liquid phenomenology.^ One can thus ex- 
pect that electron correlations will also affect the perfor- 
mance of spintronic devices made of ID wires i^iiSiiii In 
particular, local inhomogeneities, which are important to 
achieve controllable spin polarization in the Streda and 
Seba setup, strongly suppress the conductance in Lut- 
tinger liquidsji^ii^ It is thus mandatory to investigate 
how the spin polarization is affected by Luttinger liquid 
physics. This is the second goal of our work. A first step 
in analyzing the performance of the Streda and Seba spin 
filter in the presence of the electron-electron interaction 
was taken in Ref. [3- 

Standard methods such as the self-consistent Hartree- 
Fock approximation do not capture the Luttinger liquid 
physics of our setup and are known to lead to severe arti- 
facts if being applied to low-dimensional electron systems 
with Coulomb interaction. We therefore use an approxi- 
mation which is based on the functional renormalization 
group (fRG) approach to treat the two-particle interac- 
tion in a model Hamiltonian. In the absence of SOI, it 
was shown to be a reliable tool to calculate the linear 
conductance of inhomogeneous quantum wires for weak 
to intermediate interact ions j ■ "■^ ^ 

We show that for vanishing two-particle interaction 
spin polarized currents can be obtained using a similar 
mechanism as in Ref. 0, even in the presence of nonmag- 
netic leads with vanishing SOI and zero magnetic field. 
We then include the Coulomb interaction in our analysis. 
We first consider short wires (several tens of nanometers) 
in which Luttinger liquid physics does not become appar- 
ent and investigate how the spin polarization is affected 
by the two-particle interaction. We find that it strongly 
modifies the energy regime (energy of the incoming par- 
ticles) in which spin polarization can be achieved. Next, 
the focus is on system sizes for which Luttinger liquid be- 
havior is apparent in the absence of SOI (wires of several 
hundreds of nanometers). In particular, we analyze if the 
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spin polarization of an inhomogeneous wire is suppressed 
as a function of an infrared energy scale. We find that 
although the total conductance shows a power-law sup- 
pression in the presence of SOI (similar to the situation 
with vanishing SOI) the polarization does not follow such 
a scaling law. We present indications that the degree of 
spin polarization in the presence of Coulomb interaction 
might even exceed the one obtained for vanishing two- 
particle interaction. 

While in most studies on spintronic devices the corre- 
lations are neglected even if the suggested setups contain 
ID quantum wires our results clearly reveal the impor- 
tance of the two-particle Coulomb interaction in the spin 
filter suggested in Ref. @- We here refrain from making 
direct contact to existing or future experiments as we are 
mainly interested in studying the basic physics within a 
simplified model. It must certainly be extended to be 
considered realistic. However, our parameters are taken 
from a physically sensible range (see below). 

This paper is organized as follows. In the next sec- 
tion, we introduce our setup and lattice model. The tech- 
niques to obtain the linear conductance are described in 
Sect. mil In Sect. lIVl we first present our results for spin- 
polarized transport at vanishing two-particle interaction 
in the presence of leads. We next include the Coulomb 
interaction and study its interplay with the SOI, the 
magnetic field, and external potentials. Our results are 
summed up in Sect. |Vl 

II. THE MODEL 
A. Spin-orbit interaction 

The two prototypical experimental systems for strictly 
ID quantum transport are carbon nanotubes and con- 
fined electron gases which form at the interface of prop- 
erly designed semiconductor heterostructures. In the lat- 
ter, the confining potential generically leads to a siz- 
able SOI and these are the type of systems we have in 
mind in the following. But also in the former the SOI 
seems to be surprisingly large.— We choose our coordi- 
nate system such that the heterostructure confines the 
two-dimensional electron gas in z-direction and the ex- 
ternal potential results in a confinement of the electrons 
in y-direction. Therefore, the electrons are able to move 
in x-direction only. Since we assume the confinement to 
be very sharp, the different electronic subbands will be 
well separated. For a sufficiently low electron density, we 
can thus focus on the lowest subband and neglect any 
subband mixing. 

The sharp confinement leads to large electric fields, 
which induce a spin-orbit couplingi^ 

ffso = --4^<T- [ex (p--a)1 (1) 

with the electric field E = —W/e (e < is the elec- 
tron charge) being the gradient of the ambient poten- 



tial. For the ID infinite noninteracting continuum model 
(ID electron gas), SOI results in a horizontal splitting of 
the quadratic electron energy dispersion e(fc, s) with wave 
number k and s = ± being an additional quantum num- 
ber. Within a certain parameter regime a magnetic field 
in a;-direction, which couples to the electron spin, leads to 
a deformation of the lower parabolic branch resulting in 
a double well form (see the low-energy region of the two 
central dispersions of Fig. [1]). Moreover, it was shown 
in Ref. and extensively discussed in Ref. [l^ that the 
spin expectation values become fc-dependent and display 
a rich behavior in the presence of both SOI and magnetic 
field. Applying an additional step-like potential and con- 
sidering electrons with energy in regions of only double 
degeneracy leads to the above mentioned spin-polarized 
transport i^ii^ 



B. Lattice system 

The starting point of our investigation is a ID non- 
interacting tight-binding lattice Hamiltonian. As our 
single-particle basis, we choose Wannier states \ j, a) with 
j G {1, A^} labeling the lattice site and a [ denot- 
ing the spin. The spin quantization is chosen along the 
z-direction. The Hamiltonian can be written as 

^^0 = -fffroc + ^^pot + Hr + Hz , (2) 

with the free part 

N-l 
j=l a 

describing the hopping of amplitude t > and the exter- 
nal potential (e.g. due to nano-device structuring) 

N 

^p°t = ^o.A.o^o^- ■ (4) 

j=\ a 

Here, c| ^ denotes the creation operator of an electron at 
site j with spin a. In a lattice model the effect of SOI due 
to confinement in z-direction is accounted for by a spin- 
flip hopping,^*' and the effect of SOI due to confinement in 
y-direction by an imaginary spin-conserving hoppingJ^ 
Thus the SOI (Rashba) term reads 

Hr = " 51 XI '^^J (4+l,<T Cj,<T' + H.c.) (5) 

N-l 

+ H 51 "y-J (4+i.<T (*^-)<T,^' Cj.a' + H.c.) . 

The SOI coupling constants a^.j > and Uyj > are 
assumed to depend on the bond (j, j -I- 1) consideredi^ 
The effect of a magnetic field is captured by the Zeeman 
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term 



N 



Hz = 7^^Bj4_^ ^(cr:r)^^^-sin6'cos(^ (6) 

j — 1 <7,a' 



+ K)a.<T' sin 61 sin 93 + {a^)^ , cost 



with a site dependent magnetic field in (for now) arbi- 
trary direction B j — Bj (sin 9 cos tp, sin 6 sin cos 0) , and 
7 being the Zeeman couphng constant. For iV — > oo 
this lattice model was shown to give a similar low-energy 
dispersion as the continuum model as well as a similar 
energy dependence of the spin expectation valuesJ^ 

We supplement our model by a local site-dependent 
Coulomb interaction Ui,j 



N 



(7) 



■j = l a,a' 



and a bond-dependent nearest-neighbor Coulomb inter- 
action U2,j 



Af-l 



H2 = 5m^2j4+i,<TCj+i^'^4<T'Cj^'-' 



(8) 



(extended Hubbard model). The total Hamiltonian of 
the quantum wire is given by 



a — Hq + Hi + H2 ■ 



(9) 



At the sites 1 and N the ID wire is (end-)coupled to 
two semi-infinite leads. Having in mind a possible exper- 
imental realization, we assume that the SOI in the leads 
is weak and can be neglected. Furthermore, the magnetic 
field is restricted to the wire region. We assume that after 
entering the leads the electrons are independent (Fermi 
liquid behavior in higher-dimensional systems) . For a lo- 
cal lead- wire coupling and in the low-energy limit, only 
the leads' density of states at the end of the leads and at 
the chemical potential /i (energy of incoming electrons) 
matters. For simplicity, we thus model the leads as semi- 
infinite ID tight-binding chains (x = L, R) 



H. 



lead 



d'^ d 



H.c. 



(10) 



with ^ being the creation operator of the leads, j = 
—00, for the left lead, and j = N -\- 1, 00 for the 
right one. To prevent a proliferation of parameters we 
assume equal leads, set — tj^ — t, and measure all 
energies in units of the lead hopping t — 1. Our energy 
unit is therefore of the order of 1 eV. At the contacts, 
the electrons can tunnel in and out of the wire and the 
Hamiltonian of the wire-lead coupling is given by 



J2 [tLc\^^do^a + iflrfjv+l^ffCAT,^ -t- H.c] (11) 



The site- and bond-dependence of a^, az, the magnetic 
field, and the interaction matrix elements allows us to 



adiabatically turn on and off these couplings over mi 
lattice sites/bonds close to the wire-lead contacts. This 
will be done in order to prevent any unwanted electron 
backscattering from the contacts and is reminiscent of the 
gradual confinement to the ID geometry in heterostruc- 
tures. We emphasize that the precise shape of the weight 
function with which the couplings are turned on and off 
does not have any significant effect on the results as long 
as it is sufficiently smooth. In the bulk of the wire, these 
parameters reach constant values (for details, see below). 



III. METHODS 



The linear conductance 



Using the Landauer-Biittiker approach, — one can ex- 
press the spin-dependent linear conductance Ga-,^' for 
vanishing two-particle interaction in terms of the trans- 
mission T^o-' (e) 



Ga 



(^)r/'(e)de 



(12) 



with /(e) being the Fermi function. The indices a, a' 
denote the z-component of the electron spin before en- 
tering and after leaving the quantum wire, respectively. 
The spin is conserved outside the quantum wire, since 
we neglect any SOI in the leads as well as spin relax- 
ation. Using the Feshbach projection;^ it is easy to show 
that the transmission is connected to the (1,-/V) matrix 
element of the retarded single-particle Green's function 
Q'^'"' (e + iO) of the entire system (including the leads) 



2tLtR sin(fce)t/5^;^ (e - 



iO) 



(13) 



with fcg — arccos(— e/2). 

At T = 0, on which we focus from now on, the deriva- 
tive of the Fermi function is a (5-function and Eq. (fT2l) 
simplifies to 



Gfj^a' 



(14) 



with the chemical potential given by the Fermi momen- 
tum, = — 2cosfci?. For noninteracting leads, this rela- 
tion holds even if the two-particle interaction in the wire 
is finite^ where Q in Eq. (jl3p is the interacting Green's 
function of the entire system. To compute the latter we 
first integrate out the leads by projection (see the next 
subsection) and treat the remaining interacting system 
of size N using the approximate fRG procedure. 

Using standard Feshbach projection, the effect of the 
leads can be cast in energy dependent contributions to 
the (1,1) and (A^, A'^) matrix elements of the self-energyi^ 
They read 



(Slcad)i;r (z) = 5a,a'tig{z) 
(Sload)^°^r (2;) = ^<7.a'tR9(z) 



(15) 
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with the Green's function 

g{z) = z + ^T v/(2 + m)^ -4 (16) 

of the semi-infinite leads taken at the first lattice site. 
The sign must be chosen such that lim ga.a-iz) — 0. In 

\z\ — ^oo 

the following, we absorb the contribution of the leads to 
the self-energy into the noninteracting propagator Qq of 
the wire. The full interacting Green's function is given 
by the Dyson equation 

g = {g^'-^y\ (17) 

B. Functional renormalization group 

We briefly describe the fRG method used here-^ to 
approximately compute Q^i21 Detailed accounts of the 
application of this method to inhomogeneous, interacting 
quantum wires have been given in the last few yearst^^^ 
Including the SOI and the magnetic field does only re- 
quire minor extensions. 

Starting point of the fRG scheme is the propagator 
Qq which follows from the noninteracting Hamiltonian 
Eq. (21) and the self-energy contribution of the leads 
Eq. It is supplemented with a cutoff A such that 

all modes with Matsubara frequencies |w| < A are sup- 
pressed, i.e. 

g^{ico)^ei\Lo\~A)goiiLo)., (is) 

where A runs from cx) to 0. Inserting Qq in the generating 
functional of the one-particle irreducible vertex functions, 
one obtains an infinite hierarchy of coupled differential 
equations for the vertex functions by differentiating the 
generating functional with respect to A and expanding it 
in powers of the external fields. To obtain a manageable 
set of flow equations, this hierarchy must be truncated. 
In a first step, we neglect the three-particle vertex F^, 
since it is zero at A = oo and is generated only from 
terms of third order in the two-particle vertex F^, which 
are small as long as F^ does not become too large. 

For arbitrary local Ui, nearest- neighbor interaction 1/2, 
and filling n G (0, 2) of the band, the flow of the two- 
particle vertex must be kept at least to lowest (that is 
second) order to correctly obtain the scaling behavior (for 
large N) of correlation functions to leading order in the 
interaction. In the absence of SOI this is already known 
from the so-called "g-ology" modeli^ An RG analysis of 
this model shows that the two-particle scattering with 
momentum transfer 2kp of electrons with opposite spin, 
the so-called gi^± term, flows to zero. If one is interested 
in correlation functions of large systems, this scaling must 
be captured by any sensible approximation. For the ex- 
tended Hubbard model and vanishing SOI, this has been 
done in Ref. TgI using fRG. As the parameterization of the 
two-particle vertex used there relies on spin conservation, 
it cannot easily be extended to the present situation with 
SOI. Here, we are thus forced to proceed differently. 



In the extended Hubbard model the coupling gi^± is 
given by 

gi,± = Ui + 2U2 COS {2kF) . (19) 

If it is zero initially, it will not get generated in a lowest- 
order RG analysis of the corresponding "g-ology" model 
(for vanishing SOI) . If one is interested in the behavior of 
correlation functions to leading order in the interaction, 
the flow of the vertex can then be neglected altogether. 
A vanishing gi^± is achieved, if we stick to parameters 
Ui, U2, and kp such that the right hand side of Eq. (fig]) 
is zero. Neglecting terms of order U^a, where U stands 
for either Ui or U2 and a for either or ay the same 
reasoning holds if the SOI is included. When studying 
the effect of the two-particle interaction for large systems 
(Luttinger liquid behavior in the scaling limit), we thus 
exclusively consider the case 

C/2 = -C/i/[2cos(2fcF)] (20) 

and neglect the flow of the two-particle vertex. This as- 
sumption does not affect the asymptotic power-law scal- 
ing of the conductance for vanishing SOI.^'' We expect 
the same to hold for nonzero SOI and thus focussing on 
this situation does not present a severe constraint for our 
purposes. An fRG analysis of the flow of the two-particle 
vertex in the presence of SOI which would allow to in- 
vestigate novel phases resulting from the interplay of the 
two-particle interaction and the SOI^^ is left for further 
studies. 

For wires of the order of a hundred lattice sites, the 
flow of the components of the two-particle vertex is cut 
off on fairly large energy scales (oc 1/N), affects cor- 
relation functions (such as the single-particle Green's 
function we are aiming at) only weakly, and can thus 
be neglected even if one chooses parameters such that 
Eq. (j20p does not hold. Our approximation contains at 
least all terms of first order perturbation theory in the 
two-particle interaction. ^^'^^ 

Within these approximations the self-energy (neg- 
ative of the one-particle vertex) becomes frequency 
independen t -'^^1 , and its flow equation reads 

^ Ee'"°"^2VMri'.2';l,2, (21) 
uj=±A 2,2' 

where the indices 1, l',2,2' label the quantum numbers 
j, a and 

e^(zc^)= [^o-i-S^]"' . (22) 

One starts the numerical integration at a large initial 
value Aq ~ 10* and integrates down to A = 0. Following 
the description in Ref . llSl . we add a one-particle potential 
ly to the Hamiltonian Hi and H2, such that the starting 
value of the self-energy, which accounts for the finite con- 
tribution resulting from the integration of Eq. (1211) over 
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the interval (cxd,Ao], vanish. The initial condition for 
Ag ^ oo then reads 



= 



(23) 



The self-energy at the end of the flow can be con- 

sidered as an approximation to the full self-energy. Using 
the Dyson equation (flTll . the Green's function entering 
the expression for the T = linear conductance can be 
computed. 

For our model with only local and nearest neighbor in- 
teractions the self-energy is a tridiagonal matrix in the 
Wannier basis states, with each entry being a 2 x 2 matrix 
in spin space. As E is frequency independent within our 
approximation, its matrix elements can be interpreted as 
the interaction induced renormalizations of the parame- 
ters of the noninteracting model. Depending on the ma- 
trix element considered these are the SOI and spin con- 
serving hoppings, magnetic fields as well as on-site poten- 
tials all being position dependent. This implies that 
to approximately obtain the T = conductance of an in- 
teracting system the Green's function of an effective non- 
interacting problem with renormalized parameters must 
be computed. This is similar to the Hartree-Fock approx- 
imation. However, we emphasize that in our approximate 
fRG a different class of diagrams is resummed, this way 
avoiding the artifacts known to emerge if the Hartree- 
Fock approximation is used for problems of interacting 
electrons in low dimensions. 

For half filling (n — 1), the electron correlations drive 
our system towards a Mott insulator state (details de- 
pend on the parameters Ui and U2), which becomes of 
relevance for sufficiently large N.^ As we are here inter- 
ested in spin dependent transport through metallic wires 
we only consider fillings away from n = 1 when studying 
large TVi^ 



IV. RESULTS 

Following the steps described in the last section, we 
calculate the spin resolved linear conductance Go-.o-' with 
CT, <T = t, I (with respect to the ^-direction) , as a function 
of the system's chemical potential /i (energy of incoming 
electrons) and the system size N. The total conductance 
is given by the sum of the four components 

Gtotal = G|| + G|x + Gxi + Gil ■ (24) 

As our leads are free of SOI, the definition of the spin 
polarization does not involve the spin expectation value. 
This has to be contrasted to the definition of Ref. for 
a system without leads. Due to the choice of the z-axis 
as the spin quantization axis, the spin polarization in z- 
direction can be defined as the "normalized" difference 
between the probability that an electron enters the right 
lead with spin up and the probability that it enters with 
spin down 

Gfi + Gil ~ ~ ^11 




left lead | quantum wire (SOI, V and B ) | right lead 



FIG. 1: (Color online) Sketch of the system with SOI and 
parallel magnetic field in the wire. The leads have a cosine- 
like dispersion, whereas the "local" dispersion in the quan- 
tum wire has two nondegenerate branches [s = -I- (solid line), 
s = — (dashed line)]. A potential step which is turned on 
smoothly at the left contact and turned of sharply in the mid- 
dle of the wire is indicated by the dashed-dotted line. The 
energies of the incoming and outgoing electrons are confined 
within the dotted lines. 



The conductance components and spin polarization in 
the X- and y-direction can be obtained by a simple base 
transformation. The transmissions which need to be in- 
serted into Eq. p4|) are 





= C^TT 




+ ^iT 


+ ^u)/2, 


n-ix) 

■'n 


= ("^TT 


- 


+ ^iT 


-^u)/2. 




= (^TT 




-^iT 


-^ii)/2, 


q-ix) 

^11 


= {Til 






+ ^u)/2 



and 



iv) 



TT 



(y) 



n 



(y) 



it 



(y) 



U 



(26) 

{Tii~iTii+tTii+Tii)/2, 
{~iTii+Tii+Tii+iTii)/2, 
{iTii + Til + - i^ii) /2 , 
{Tii+iTii-iTii+Tii)l2, (27) 



P. 



G 



total 



(25) 



where T„^a' denotes the transmission with respect to the 
z-direction. The corresponding polarizations follow as in 
Eq. (ESI) with G replaced by G^"^) and G'^y\ respectively. 



A. Vanishing Coulomb interaction 

1. A potential step 

To investigate whether spin polarized currents can be 
achieved using a mechanism similar to that of Ref. in 
the presence of nonmagnetic leads, we consider a setup 
that is close in spirit to the one studied there. We there- 
fore add a potential step to the wire which can be used to 
shift the energy into the region of only doubly-degenerate 
states as sketched in Fig. [TJ To prevent any backscatter- 
ing at the left contact it is turned on smoothly there 
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FIG. 2: (Color online) Spin polarization in x-, y- and z- 
direction as a function of fi for a homogeneous system [Vo ~ 0, 
a)] and a system with potential step [Vo = 0.4, b)] for sys- 
tem parameters A'^ = 100, mi = 20, II — ta = l,ay = 
0, Qz = 0.8, and yB^ = 0.6. Inside the region with two 
conducting channels, there is no significant spin polarization. 
The high and low energy polarization is independent from 
Vo in good approximation. The potential leads to effects for 
M e [-2 + -2 + jB^ + Vo]. 
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FIG. 3: (Color online) Conductance and spin polarization for 
a system of N = 101 lattice sites (mi — 20) with a single 
impurity Vi in the middle of the system in presence of a mag- 
netic field yBx = 0.6 and different SOI couplings {ay = 0). 
a) The homogeneous system shows a perfect step-like polar- 
ization in a;-direction. b) The shape of Px is washed out in 
presence of a strong single impurity Vi = 4 and the conduc- 
tance is severely decreased, c) Small SOI coupling az = 0.2 
leads to finite Py, Pz and polarization oscillations, d) The 
oscillations are enhanced by increasing a^. 



and turned off sharply in the middle of the wire (dashed- 
dotted line in Fig.[T]). Due to the smooth variation of the 
system parameters at the left contact, the cosine-like dis- 
persion of the leads is "adiabatically" transformed into 
the "local" dispersion in the wire. A similar "adiabatic" 
transition occurs at the right contact. As the achieved 
wire dispersion is an essential ingredient of the spin filter 
to perform properly^^ii^ modelling a smooth variation at 
the contacts is mandatory. However, assuming a gradual 
crossover from higher-dimensional leads to the ID wire 
appears to be quite natural in heterostructures. 

Fig. [2] shows the spin polarization Pj, i = x,y, z, for a 
homogeneous system [Vq = 0, Fig. [2^)] and a system with 
potential step [Vq = 0.4, Fig. [^b)] for system parameters 
TV = 100, mi = 20, ay = 0, = 0.8, and -/B^ = 0.6. 
We chose the value for Uz to be larger than the spin- 
orbit parameter a ~ 0.1 eV extracted from Ref. |30| for 
bulk InAs, since the SOI is significantly increased in 
semiconductor heterostructures. The lead- wire tunnel 
contacts modelled by ii/zj are assumed to be "perfect" 
= tj^ = 1, a situation on which we focus from now 
on. The curves show a large spin polarization in x- and 
y-direction for fi ^ [—2 + jB^ + Vo,2 — 'fB^], i.e. as long 
as there is only one conducting channel. This polariza- 
tion is a "trivial" band effect due to the Zeeman splitting. 
Without SOI, one obtains a perfect spin-polarization in 
cc-direction, — ±1, in this interval [see Fig. [3^)]. For 
finite SOI, the spin is rotated out of the a:-direction lead- 
ing to non- vanishing components Py and Pz- Similarly 
to the continuum situation^i^ {(^y) niainly leads to a 
spin rotation into the y- (z-) direction. 

The most interesting energy regime in connection with 
Ref. ^is ii&[-2 + -fB^, -2 + jB^ + Vq] in Fig. Wp)- In 
this interval of width Vq, we observe nonvanishing, but 



strongly oscillating components of the polarization. The 
oscillations can be traced back to those of the conduc- 
tance components and are absent in the continuum model 
of Ref. [y. By a mechanism similar to the one discussed 
in Ref. 0, we obtain a nonvanishing polarization. As was 
pointed out in Ref. [l^ for the lead- less continuum model, 
the polarization P^ (in x-direction) only depends on the 
absolute value a and not on the direction of the effective 
SOI field and \Py/Pz\ = \az/ay\. Averaging over the 
polarization oscillations, this holds in very good approx- 
imation also for the lattice model with leads, although 
the polarization here is defined by Eq. (j25p and not in 
terms of spin expectation values Because of the large 
oscillations (absent in the lead-less continuum model) the 
spin polarization reacts very sensitive to changes in the 
chemical potential /i. This could be of interest to control 
the spin polarization in future spintronic devices. 



2. A single impurity 

Another interesting polarization effect occurs, if we in- 
sert a localized impurity of strength Vi in our system 
with SOI and parallel magnetic field. The exact position 
of the impurity within the bulk part of the wire does not 
matter and we here locate it in the middle of the wire. 

Fig. [3] shows the conductance and spin polarization of 
a system with N — 101 lattice sites (toi ~ 20) with a 
single impurity Vi for a magnetic field jB^ = 0.6 and 
different SOI couplings az {ay = 0). The system with- 
out SOI and without any impurity {Vi — 0) shows a per- 
fect step- like polarization in x-direction ("trivial" band 
effect), whereas Py = Pz = [see Fig. [3^)]. The step 
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is smeared out and Px can be tuned smoothly in pres- 
ence of an intermediate to strong single impurity Vi = A 
while the total conductance is severely decreased in this 
case [see Fig. For small SOI, = 0.2, P^, is stih 

the dominant polarization component, but Py and P^ 
become finite. All polarization components show oscilla- 
tions for ^ S [—2 -|- jBx, 2 — jBx], which become heavily 
pronounced at the edges of this interval [see Fig. ^)]. 
For large SOI [see Fig. [3ji)], — 0.8, we observe the 
same behavior for /i ^ [—2 -I- — ^Bx\ as in the 

impurity-free case [see Fig. [2^)] with P^ playing only a 
minor role. The oscillations of the polarization compo- 
nents are more pronounced compared to the case with 
small ckz, especially for P^- Moreover, the total conduc- 
tance is enlarged due to the smaller ratio of Vi and the 
effective hopping yjt"^ + al- Due to the large oscillations, 
each spin polarization component reacts very sensitive to 
changes in /i and can therefore be tuned by adjusting /i. 

In this section we studied systems with TV ^ O(IO^) 
lattice sites, corresponding to wires of the order of tens 
of nanometers. For vanishing two-particle interaction in- 
creasing the number of lattice sites does not affect the 
results obtained here. 
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FIG. 4; (Color online) Total conductance and spin polariza- 
tion as a function of the chemical potential ^ for a system 
with A'^ = 100 lattice sites (mi = 20), potential step Vo = 0.4, 
magnetic field 'yBx =0.6 and different Coulomb interaction 
U = U\ = 2U2- Increasing the Coulomb interaction [from a) 
to b)] for constant SOI Qz ~ 0.8 {oy — 0) leads to a broaden- 
ing of the energy regime with two open conducting channels. 
The spin polarization vanishes in this regime, but reveals a 
strong dependence on jj, in the regime with only one conduct- 
ing channel. 



B. The effect of the Coulomb interaction 

We now add the terms Hi and H2 given in Eqs. ([7]) and 
([5]) to our Hamiltonian and use the fRG to approximately 
compute the conductance. We first study short quantum 
wires with N ~ 0(10^) lattice sites and investigate how 
the energy regime (of the incoming electrons) in which 
spin polarized currents can be obtained is modified by the 
Coulomb interaction. Changes can be traced back to the 
interaction induced renormalization of the parameters of 
the noninteracting model. '^^ 

In a second step we focus on a fixed chemical potential, 
at which spin polarization is observed and study how this 
is modified if the system size N is increased. The inverse 
of N, more precisely S oc vp/N (with the Fermi velocity 
vp), presents an infrared energy scale in our setup. In 
the absence of SOI and a magnetic field, inhomogeneities, 
such as single impurities and potential steps, are known 
to lead to a power-law suppression of the conductance as 
a function of infrared energy scales J^ii^ii^ii^ We demon- 
strate that while the total conductance shows a similar 
behavior for finite SOI, the polarization does not display 
scaling behavior. 

To model a gradual transition from a higher- 
dimensional to a ID system at the lead-wire contacts, we 
gradually increase Ui and U2 over mi lattice sites starting 
at the contacts. This prevents electron backscattering at 
the contacts due to the inhomogeneous two-particle in- 
teraction and one achieves unitary total conductance in 
the absence of an external single-particle potential. As 
the details of the variation of the interaction do not mat- 
ter as long as it is sufficiently smooth,— we chose the 
same weight function as for the single-particle parame- 



ters. For an increasing chain length, the number of lat- 
tice sites over which the two-particle interaction and the 
single-particle parameters are turned on and off close to 
the contacts must be increasedi^ 



1. Short wires 

Fig. m shows the total conductance and spin polariza- 
tion as a function of the chemical potential fi for a system 
with the same parameters as in Fig. [^Jd) but for nonvan- 
ishing Coulomb interaction. We consider a constant ra- 
tio of local and nearest-neighbor interaction, U1/U2 = 2 
(which appears to be a rather physical value), denote the 
local interaction Ui by U and study its affect on Gtotai 
and Pi, i — x,y,z. Comparing Figs. Eb),!!^) andHb), 
one notes that an increase of U leads to a broadening of 
the energy regime with two open conducting channels, 
that is "perfect" conductance Gtotai ~ 2e^//i, and van- 
ishing spin polarization. The Coulomb interaction thus 
leads to a decrease of the energy range in which spin 
polarization can be achieved. We emphasize, that within 
our approximation the above behavior can be understood 
in terms of transport through a quantum wire with van- 
ishing two-particle interaction, but renormalized single- 
particle parameters (regular and spin-flip hoppings, mag- 
netic field, and on-site potentials all of them depending 
on the lattice site index j) given by the self-energy at the 
end of the fRG flow^ 

We now investigate the influence of the Coulomb in- 
teraction on the transport properties of a short wire 
with a single impurity Vi in the middle of the quan- 
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FIG. 5: (Color online) Total conductance and spin polariza- 
tion components for a system with A'^ — 101 lattice sites 
(mi = 20), SOI coupling parameters ay = Oz = 0.5, mag- 
netic field in i-direction, ^Bx = 0.6, and a single impurity 
Vi = 2 in the middle of the system. Increasing [from a) to 
b)] the two-particle interaction U leads to an increase of the 
energy regime with two conduction channels, but an overall 
decrease of the total conductance. 



turn wire. Fig. [5] shows the total conductance Gtotai 
and polarization components Pi, i = x,y,z, for a sys- 
tem with N = 101 lattice sites (mi = 20), SOI param- 
eters Uy — ^ 0.5 and magnetic field 'yBx — 0.6 in 
x-direction for Vi = 2. Increasing the Coulomb inter- 
action, we observe two tendencies. The first one is the 
increase of the energy regime with two conduction chan- 
nels already found in the presence of the potential step. 
As the electrons scatter off the single impurity this does 
not lead to unitary conductance in the present case (in 
contrast to the case of a potential step). The second ef- 
fect is an overall decrease of the total conductance with 
increasing U . This was to be expected as it is known that 
the effective strength of inhomogeneities increases in the 
presence of Coulomb correlations, eventually leading to 
power-law suppression for sufficiently large N (that is 
sufficiently small energies J) , 12,13,15^16 

For the fairly short wires studied in this subsection, 
results from first order perturbation theory in the two- 
particle interaction (for the self-energy) lead to qualita- 
tively similar results as those obtained by our approxi- 
mate fRG scheme. This does no longer hold for the longer 
wires studied next (Luttinger liquid behavior). 



2. Luttinger liquid behavior in long wires 

In the absence of SOI and a magnetic field, the lin- 
ear conductance of a Luttinger liquid wire with a single 
impurity Vi is known to show power-law suppression as 
a function of an infrared energy scale (e.g. the temper- 
ature or (5 oc vf/N), provided the latter is sufficiently 
small (scaling regime) J^ii^ii^ii^ Surprisingly Gtotai van- 
ishes in the asymptotic low-energy limit even for small Vi . 
However, the energy scale beyond which scaling holds de- 



FIG. 6: (Color online) Dependence of the scaling exponent on 
the strength of the SOI (varying Oy as well as Oz) for U — 0.3 
and filling n — 0.75. The circles are the data extracted from 
fitting Gtotai (5) by a power law for systems of a few thousand 
lattice sites. The error is of the order of the symbol size. 
The effective hopping tcs depends on ay as well as and is 
defined in Eq. (|30p . The line shows the analytic expression 
Eq. (|28p with t t^s- The systematic deviation of the fRG 
data from the analytic expression can be explained by higher 
order corrections in U included in the numerical data, but not 
in Eq. 



pends on Vi and becomes exponentially small for small 
Vi , implying that exponentially long chains must be stud- 
ied to observe power-law behavior. The scaling exponent 
is a function of the interaction and the filling, but is inde- 
pendent of the strength of the bare impurityJ^'i^'i^ii^ For 
the extended Hubbard model and our choice of the ratio 
U1/U2 = U/U2 = —2cos{2kp) (vanishing two-particle 
backscattering gi±; see Sect. IIIIBp it is given by 



2a, 



/i^ — 4 cos(7rn) 
(2 - /i2)(27rtsin(7m/2)) 



U 



(28) 



to leading order in U.— It was shown earlier that the 
approximate fRG procedure captures this power-law be- 
havior and correctly reproduces the scaling exponent to 
order U In order to be able to compare our results 
to Eq. (P5|) we need to know n and therefore tune the 
additional one-particle potential viU,^), introduced in 
Sect, nil Bl such that the filling of the ID quantum wire 
with electrons in presence of Ui and U2 corresponds to 
the filling n — 2 arccos (— /j,/2)/7r of the noninteracting 
leads at given ^. Following Ref.[l^ the starting values of 
the self-energy at lattice site j due to integration of the 
flow equations from 00 down to Aq are given by 

-/^^ - (\ui,j + 2U2^ <J,,,'<5.,.,, J e {2, . . . , TV - 1} , 



- f \Ui,, + U2,, ) j = l,N , (29) 



3,3 



up to corrections of order l/Ag. 

We first investigate if power-law scaling in the pres- 
ence of a single impurity is also found for nonvanishing 



9 




0.15 



500 



N 



1000 



FIG. 7: (Color online) Main plot: Dependence of the polar- 
ization on the system size N for for Vi = 8, = = 0.5, 
yBx = 0.6, n — 0.9, and U = 0.8. Inset: The total conduc- 
tance as a function of TV on a log-log scale. The parameters 
are as in the main plot. 



SOI, but zero magnetic field, a situation in which no 
spin polarization is found. To this end we compute the 
total conductance as a function of S oc vp/N for fixed 
(small) ?7 ~ 1 (range of applicability of our approximate 
fRG procedure), fixed filling n ^ 1 (see the discussion in 
Sect, nil B|) . fixed SOI, and fixed intermediate to large Vi 
[such that the scaling regime is reached for N ~ O(IO^)]. 
The data for GtotaK^) can be fitted by a power-law and 
we extract the asymptotic exponent. An example of the 
power-law behavior (as a function oi N ^~^) in the 
case of an additional magnetic field is shown in the inset 
of Fig. [71 The scaling exponent depends on the SOI via 
an effective renormalized spin conserving hopping 



*off 



at, 



(30) 



that is, the dependence of 2as is given by Eq. ([281) with 
t replaced by tcff. This is shown in Fig. [6lfor U = 0.3, 
n = 0.75, and different ay as well as a^. The systematic 
deviation of the data from the analytic expression ([28]) 
with t —> tcff can be explained by higher order corrections 
in U included in the numerical data, but not in Eq. (j28p . 
The error of the fRG exponents extracted from fitting 
Gtotai(<5) for systems of a few thousand lattice sites is of 
the order of the symbol size. 

Next we consider the case of a single impurity in the 
presence of SOI and a magnetic field in x-direction, a 
situation with nonvanishing polarization (see the discus- 
sion in Sec. IIV A 2p . The main part of Fig. [71 shows the 
TV ~ dependence of the three components of the po- 
larization for Vi = 8, ay = az — 0.5, jB^ = 0.6, n = 0.9, 
and U = 0.8. In the inset we present the corresponding 
total conductance on a log-log-scale. While the latter 
shows a clear indication of power-law suppression, the 
polarization oscillates with a constant amplitude Aosc- 
A similar oscillation is found for [7 = 0. In fact, the 
amplitude of the oscillation shows a nonmonotonic de- 
pendence on U. Starting from U = 0, it first decreases 




0.05 



FIG. 8: (Color online) Amplitudes Aosc of the oscillations of 
the three components of the polarization (see main part of 
Fig. [7} as a function of the two-particle interaction U for the 
same parameters as in Fig. [71 Because of Oy = a^, the ampli- 
tudes of the polarization in y- and z-direction are supposed 
to be the same. The discrepancy of these two data sets thus 
gives an estimate how accurately the amplitude can be read 
off from the numerical data. 



linearly with increasing [/ up to [7 = 0.3, but increases 
for U > 0.3 (see Fig. HI). For U = 0.8 as shown in Fig. [7| 
the amplitude is roughly a factor of 1.5 larger than in the 
noninteracting case. Fig. \5\ indicates that the details of 
this behavior depend on the filling (chemical potential of 
the incoming electrons). The above result implies that 
although the total current through (total linear conduc- 
tance of) a quantum wire with a single impurity, SOI, and 
a magnetic field in the direction of the wire is generically 
strongly reduced as a function of TV in the presence of 
the Coulomb interaction, the degree of spin polarization 
of the current stays constant. We note that the increase 
of Aosc as a function of U (see Fig. [Q) has to be con- 
sidered with caution as our approximation is only valid 
for sufficiently small U, while the increased polarization 
(compared to the noninteracting one) requires a small to 
intermediate U. It would thus be important to inves- 
tigate the polarization using alternative methods. We 
verified that the exponent of the power-law suppression 
of Gtotai(<5) is independent of ^B^ and thus to leading 
order in U given by Eq. ([^ with t replaced by tcff- 

It would now be very interesting to investigate the scal- 
ing of the conductance and the dependence of the polar- 
ization on the length of the quantum wire for the setup 
with a potential stepi^. Unfortunately, it turns out that 
the relevant 2fci^^-scattering component of the potential 
step which is turned on adiabatically at the left contact 
but is turned off abruptly in the middle of the wire with 
a step height of a few ten percent of the band width is 
fairly small. We thus cannot reach the scaling regime 
for the accessible system sizes of order 10^ to lO'*. For 
larger step height the physics is dominated by "trivial" 
band effects we are not interested in. As for weak (local) 
single impurities the total conductance in the presence 
of the potential step decreases weakly with increasing TV 
(no power-law scaling at the corresponding S) while the 
components of the polarization oscillate with TV with an 
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amplitude which is constant. One can enhance the effect 
of the potential step (the size of the 2/cj--scattering com- 
ponent) by adding a (sufficiently strong) single impurity 
located in the middle of the system. For this combined 
inhomogcneity we observe exactly the same behavior as 
shown in Fig. [7] for a pure (strong) single impurity. Hav- 
ing this in mind we conjecture that a potential step with 
a sufficient large step size not limited by the finite band 
width will have the same effect on the total conductance 
(power-law suppression with length of wire) and the spin 
polarization (oscillation with an amplitude independent 
of the wire length) as a single localized impurity. 

V. SUMMARY 

We have investigated the effect of SOI, a magnetic field 
and the Coulomb interaction on the transport properties 
of a ID quantum wire attached to two semi-infinite non- 
interacting and nonmagnetic leads. Motivated by ana- 
lytical calculations for a lead-less, noninteracting system 
described by a continuum model,—— we have constructed 
a corresponding lattice model. The combined effect of 
SOI and a magnetic field led to a spin polarization which 
could be varied over a wide range in the presence of inho- 
mogeneities (single impurity, potential step) by adjusting 
the energy of the incoming electrons (chemical potential). 

Using the fRG, we were able to include the Coulomb 



interaction in our system. We distinguished the cases 
of short quantum wires (of the order of a few tens of 
nanometers) and long ones (scaling regime; a few hun- 
dred nanometers). For short wires, we showed that the 
energy regime for which spin polarization can be found 
strongly depends on the Coulomb interaction and might 
even become very small depending on the other system 
parameters. For long wires the well-known power-law 
suppression of the total conductance of an inhomoge- 
neous Luttinger liquid as a function of the system size was 
obtained with a scaling exponent which depends on the 
two-particle interaction, the filling (both as for vanishing 
SOI), but also on the strength of the SOI via an effective 
nearest-neighbor hopping. However, the spin polariza- 
tion as a function of system size exhibited oscillations 
with constant amplitude, not signaling any suppression 
on low energy-scales (that is for large system sizes). We 
found indications that the amplitude of the oscillations 
and thus the degree of spin-polarization, might even be- 
come larger than for vanishing Coulomb interaction. 

Our results indicate the importance of the two-particle 
Coulomb interaction in the spin filter suggested in Ref. ^. 
In most studies on spintronic devices these correlations 
are neglected even if the suggested setups contain ID 
quantum wires. Obtaining a deeper understanding of 
the performance of ID spin filters within more realistic 
models of interacting electrons presents a challenge for 
future theoretical studies. 
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